To find evidence for the existence of grains of millimetre 
sizes in disks around in T Tauri stars, implying grain growth. 
The Australia Telescope Compact Array (ATCA) was used 
to observe 15 southern T Tauri stars, five in the constella- 
tion Lupus and ten in Chamaeleon, at 3.3 millimetre. The 
five Lupus sources were also observed with the Submillimeter 
Array (SMA) at 1.4 millimetre. Our new data are comple- 
mented with data from the literature to determine the slopes 
of the spectral energy distributions in the millimetre regime. 
Ten sources were detected at better than 3<x with the ATCA, 
with cr ~ 1-2 mjy, and all sources that were observed with 
the SMA were detected at better than 15cr, with a ~ 4 mjy. 
Six of the sources in our sample are resolved to physical radii 
of ~ 100 AU. Assuming that the emission from such large 
f^*) \ disks is predominantly optically thin, the millimetre slope can 
< ^ • be related directly to the opacity index. For the other sources, 
the opacity indices are lower limits. Four out of six resolved 
"q ■ sources have opacity indices < 1, indicating grain growth to 
/"■n millimetre sizes and larger. The masses of the disks range from 
< 0.01 to 0.08 M Q , which is comparable to the minimum mass 
solar nebula. A tentative correlation is found between the mil- 
limetre slope and the strength and shape of the 10-//m silicate 
| feature, indicating that grain growth occurs on similar (short) 
timescales in both the inner and outer disk. 
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Abstract. Low-mass stars form with disks in which the coagulation of grains may eventually lead to the formation of planets. 
It is not known when and where grain growth occurs, as models that explain the observations are often degenerate. A way to 
break this degeneracy is to resolve the sources under study. 



1. Introduction 

Disks of dust and gas are observed around many young 
stars. According to the so-called core-accretion model 
( Safr onov & Zviaginalll969l) . planetary systems such as our 
own Solar System are formed in these circumstellar disks: the 
solid particles coagulate to form larger grains, which will grow 
to eventually form planets. The grains mainly consist of carbon 
and silicates. The silicate grains are readily observed through 
their 10- and 20-jum features. Both the change from amorphous 
to more crystalline grains, as well as the growth of grains from 
submicron sizes to sizes of several yum have been observed 
by the Infrared Space Observatory (ISO, iMalfait et atll 19981 
Ivan Boeke l et al. 2005), and more recen tly by the Spitzer Space 
Telescope (iKessler-Sila ccTet alJ l2006). Although the qualita- 
tive picture of grain growth has become much clearer over the 
last few years, several quantitative details are still under discus- 
sion. Open questions include for example the timescale over 
which grain growth occurs and how this relates to the disk's 
phys ical structure (e.g., i ts temperature and density profile). 
See lDominik et alJ ll2006l) for a recent discussion of both labo- 
ratory measurements and theoretical modelling of the aggrega- 
tion of dust in protoplanetary disks. 

A large sample of solar-mass T Tauri stars have recently 
been observed with the InfraRed Spectrograph (IRS) on board 
the Spitzer Space T elescope, in the c ontext of the "Cores to 
Disks" (c2d) legacy fevans et alJ2003l) and other programmes. 
Most of the sources in the c2d sample show 10- and 20- 
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fim amorphous silicate features jKessler-Silacci et"all E006). 
confirming the results of earlier ISO and ground-based 10- 
nm o bservations (e.g. Jvan Boekel et all2003tlPrzvgodda et alJ 
2003), and extending the observed protoplanetary disk sample 
to lower mass objects. The data indicate a large variety of sil- 
icate profiles, ranging from strongly peaked silicate bands and 
steeply rising spectral energy distributions (SEDs) to "boxy" 
silicate profiles and flat SEDs. The boxy features with low 
feature-to-c ontinuum ratios are int erpreted as grain growth to 
micron size ( Bo uwman et al.l200ll) . 

One possible explanation for the different spectra and SEDs 
is that grain growth and the s hape of the disk are related. 
IPuUemond & Dominild J2004al) identified a correspondence 
between the growth of grains in a circumstellar disk and the 
evolution of the disk from a "flaring" to a "self-shadowed" ge- 
ometry. In their models, the larger, more massive dust grains 
settle to the midplane of the disks as the grains grow, and the 
initially flared disks evolve into flatter, self-shadowed disks. To 
what extent this process is related to the age of the young star is 
still under debate; there are indications that some young stellar 
objects evolve more quickly than others. 

Numerical models show that the process of settling 
and coagulation is rapid, being well underway at distances 
of 1-30 AU from the ce ntral star in ~ 10 4 yrs (e.g., 
iNomura&Nakagawal EoOfi). Models predict that the slopes 
of the SEDs in the millimetre wavelength range will be- 
come shallower as the grains in the disk grow to millime- 
tre (mm) and subsequently centimetre (cm) and large r sizes 
JPullemond & Dominild l2004at ID' Alessio et alJ 120061) . It is 
therefore necessary to observe these sources at larger wave- 
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lengths than the infrared. Furthermore, millimetre observations 
probe the entire disk, including the cold midplane, whereas in- 
frared observations can only probe the hot surface layer of the 
inner disk. However, a shallow millimetre slope in itself is not 
enough evidence for grain growth in the disks, since an excess 
flux at long wavelengths may also be caused by a very sm all, 
optically thick disk (see, e.g.. iBeckwith & Sargentl ll991). To 
break this degeneracy it is necessary to resolve the disks to 
determine their actual sizes, since if the physical disk size is 
known, a reasonable estimate of the disk's opacity can be ob- 
tained. To resolve the disks around T Tauri stars, interferomet- 
ric observations are indispensable. 

Considerable progr ess has been m ade in this field over the 
last several years. Wil ner et alJ yOOO) resolved the inner disk 
of the classical T Tauri star TW Hya in dust emission at 7 
mm using the Very Large Array (VLA - ). Icalvet et all (|2002) 
did extensive modelling of the SED of this source and showed 
that the dust grains in the disk m ust have grown to sizes of at 
least ~ 1 cm. Wil ner et. al J feoO^ used the Australia Telescope 
Compact Array (ATCA) to stu dy TW Hya at 3.4 mm and f ound 
that a passive two-layer disk llChiang & Goldreichll 19991) pro- 
vides a reasonable model to explain the observations. Ten T 
Tauri stars in the Taurus-Auriga st ar-forming regions wer e re- 
solved at 7 mm with the VLA bv iRodmann et alJ J2006h : the 
majority of these show strong evidence for grain growth to at 
least millimetre-sized dust. 

The more massive counterparts of the T Ta uri stars, the 
Herbig Ae/Be stars, have been studied extensively. M eeus et"al] 
(2001) classified 14 isolated Herbig Ae/Be stars into two 
groups based on the shape of the SED. Disks with a flared outer 
part of the disk show a rising mid-infr ared (20-100 u rn) excess 
(Group I in the classification by iMeeus et alJl200ll) . whereas 
self-shado wed disks have mo re modest mid-infrared excesses 
(Group /D. lAckeet ail fe004 compared the millimetre slopes 
for a sample of 26 Herbig Ae/Be stars, based on single-dish 
data, and found that Group II sources have in general a shal- 
lower slope, consistent with grain growth to larger sizes than in 
the Group I source s, and observationally confirmi n g the m odel 
predic tions from iDullemond & Dominikl j2004al) . lNatta et aD 
(2004) analysed interferometric observations of nine pre-main- 
sequence stars, mostly Herbig Ae stars, and find that the obser- 
vations are well explained with dust size distributions contain- 
ing boulders of up to metre sizes. 

We have used the ATCA to observe 3.3-mm continuum 
emission from a sample of 15 southern T Tauri stars, five in 
Lupus and ten in Chamaeleon. Compared with the Taurus star- 
forming cloud, the T Tauri stars in Lupus and Chamaeleon 
are generally somewhat older, and the Lupus clouds are 
muc h richer in ver y low mass stars than the other clouds 
(Hug hes etai1ll994l) . The sources IM Lup and WW Cha were 
also observed in spectral-line mode to search for HCO + as 
a tracer of the molecular gas. The observations of the Lupus 
sources were followed up with 1 .4-mm observations using the 
SubMillimeter Array (SMA). The sample and the observations 
are described in Sect. [2] The basic results are presented in 
Sect. [5] and further discussed in Sect. |4] We will summarize 
our results and draw some conclusions in Sect.|5] 



2. Observations 

2. 1. Source selection 

The ATCA was used to observe 15 southern T Tauri stars 
(listed in Table Q at 3.3 mm. The sources were selected 
to overlap with the sample observed in the c2d programme. 
We furthermore selected sources with strong 1.3-mm fluxes 
( iHenning et aljfl993t iNiirnberger et alJll997l) to improve the 
chances of detection at 3.3 mm. 

The distances used in this work are 150 + 20 pc to Lupus 
I (HT Lup, GW Lup) and L upus II (IM Lup RU Lup), 200 + 
20 pc to Lupus III (HK Lup) JComer6rj2006l) . and 160 ± 15 pc 
to Chamaeleon I JWhittetetalJll997l) . The distances to most 
of these p re-main-sequence stars are not very well constrained 
(see, e.g., IComer6iJl 2006: van Kempen e tal1l2 006. for discus- 
sions on the distances to the Lupus clouds). The distance to 
the isolated source T Cha, however, is known t o be closer at 
66^2 P c (Hipparcos. lvan den Ancker et alJl9 98). 

2.2. ATCA observations 

ATCA 1 observations were carried out in July 2003, October 
2004, and August 2005. All sources were observed at 3.3-mm 
continuum in double sideband. The primary beam of the ATCA 
antennae is * 35" at 3 mm. Each sideband consisted of 32 
channels and had an effective total bandwidth of 128 MHz. 
Furthermore, the sources WW Cha and IM Lup were observed 
in a dual mode: the lower sideband was a 512-channel band 
with an effective total bandwidth of 8 MHz to provide high- 
frequency resolution for detection of the HCO + J = 1-0 line, 
and the upper sideband was used as a wideband channel with 
again 32 channels and a total bandwidth of 128 MHz. These 
two sources were observed for a complete track to reach a sim- 
ilar RMS noise in the continuum as that for the other sources, 
and to maximize the possibility of detecting the HCO + J = 
1-0 line. The resulting velocity resolution of the narrow-band 
observations was 0.11 km s _1 , the velocity coverage was ~24 
km s . The ATCA was in the EW214 configuration at the time 
of the observations in 2003 (three antennae equipped with 3- 
mm receivers, baselines of 31-107 metres) and in the H214C 
configuration in 2004 and 2005 (five antennae, baselines of 82- 
247 metres). The data were calibra ted and reduced with the 
MIRIAD package JSault et alJll995T) . The quasars PKS 1622- 
297 and PKS 1057-797 served as gain calibrators for the Lupus 
and Chamaeleon sources, respectively, and the absolute fluxes 
were calibrated on Mars or Uranus. The calibration is estimated 
to have an uncertainty of ~20 %. The passbands were calibrated 
on the quasars PKS 0537-441 and PKS 1253-055. 

The phase centre was offset from the source position by 
5" to avoid possible artefacts. There was some overlap in the 
samples that were observed in the three runs. The 2004 data 
suffered badly from unstable weather and were not used in the 
analysis, whereas the 2003 and 2005 data were consistent. In 
2003, there were only three ATCA antennae fitted with 3 -mm 

1 The Australia Telescope Compact Array is part of the Australia 
Telescope which is funded by the Commonwealth of Australia for op- 
eration as a National Facility managed by CSIRO 
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Table 1. Source list of sources observed with the ATCA. 



Source 


Cloud 


Age 11 


Spectral 


Luminosity c 


Distance 11 


When observed c 


Other names' 






(Myr) 


type b 


CU>) 


(pc) 






SYCha 


Cha I 


3-5 


M0 


0.35 


160+ 15 


2005 


IRAS 10557-7655, SZ 3, 
















HBC 565, Cha T4 


CR Cha 


Cha I 


1-2 


K0,K2 


2.8 


160+ 15 


2003+ ; 2005 


IRAS 10578-7645, SZ 6, 
















HBC 244, Cha T8 


CS Cha 


Cha I 


2-3 


K4 


1.3 


160+ 15 


2003+ ; 2005 


IRAS 11011-7717, SZ 9, 
















HBC 569, Cha Til 


DI Cha 


Cha I 


3^t 


G1,G2 


8.9 


160+ 15 


2005 


IRAS 11059-7721, SZ 19, 
















HBC 245, Cha T26, KG 6 


KG 28 


Cha I 




K7 


1.1 


160+ 15 


2005 


SZ 22, Cha T29 


Glass I s 


Cha I 


2-3 


K4 


1.3 


160 ± 15 


2005 


IRAS 11068-7717, ChaT33, KG 39 


KG 49 


Cha I 






15 


160 ± 15 


2005 


IRAS 11072-7727, ChaIRN 


WW Cha 


Cha I 


0.4-0.8 


K5 


2.2 


160 ± 15 


2003+ ; 2005 


IRAS 11083-7618, SZ 34, 
















HBC 580, Cha T44 


XX Cha 


Cha I 


10-40 


Ml 


0.10 


160 ± 15 


2005 


IRAS 11101-7603, SZ39, 
















HBC 586, Cha T49 


TCha 


Isolated 


> 12.5 


G2 


1.3 


66!" 


2003+ ; 2005 


IRAS 11547-7904, HBC 591 


HT Lup h 


Lup I 


0.4-0.8 


K2 


6.0 


150 ± 20 


2003+ ; 2005 


IRAS 15420-3408, SZ 68, HBC 248 


GWLup 


Lup I 


1.3-3.2 


M2 


0.23 


150 ± 20 


2003+ ; 2005 


IRAS 15435-3421, SZ 71, HBC 249 


IMLup 


LupH 


0.09-0.6 


M0 


1.3 


150 ± 20 


2003+ ; 2005 


IRAS 15528-3747, SZ 82, HBC 605 


RULup 


Lup II 


0.04-0.5 


K7-M0 


2.1 


150 ± 20 


2003+ ; 2005 


IRAS 15534-3740, SZ 83, HBC 251 


HKLup 


Lup III 


0.7-1.4 


M0 


0.62 


200 ± 20 


2003 


IRAS 16050-3857, SZ 98, HBC 616 



a Ages adopted from lHughes et alJ |l994l).lLawson et al Jll996l). an d van den Ancker et all ll998f). 

b Spectral types adopted fro m |Her^i^^Kameswara Rao] Jl972h . Gauvin & Stroml i!992h. hughes et all < 19941) . lAlcalaetalJ Jl99l . 

Ivan den Ancker et all ll998{). andlComeron et alJ Il999l). 
c Luminosities adopted fr omlHughes et alJ<199 4.lLawson et alJjl9 9d).IChen et"ai]jl997l). a nd lvan den Ancker et"al]jl998l) . 
d Distances adopted from Whi ttet et alJ Jl997l) . lvan den Ancker et alJ ll998h . lComer6nl 120061) . 
c The observations marked with a t were not used in the anal ysis. 

f Catalogue names: KG = [KG2001] iKenvon & Gomej j2001[). IRAS = Infrared Astronomical Sat ellite, SZ = Southern Z wicky 

iRodeers et alll978l> . HBC = Herbig+Bell Catalog iHerbig & Belll988h . Cha T = Assoc Chamaeleon T <Henize & Mendozall973l) . 
g This source is a binary with a separation of 2'.'4; the spectrum quoted is that of component A. 

h This source is a binary in 2MASS K-band images with a separation < 3". The spectrum quoted includes both sources. 



receivers, resulting in a large and elongated beam. Five anten- 
nae with 3-mm receivers were available in 2005, greatly im- 
proving the resolution and beam shape of our observations. 
With the exception of HK Lup, which was only observed in 
2003, we decided to use only the 2005 data for our analysis. 
For an overview of the ATCA observations, see Table|2] 

2.3. SMA observations 

The five Lupus sources in our sample were observed in double- 
sideband continuu m (217 and 22 7 GHz, primary beam ~ 55") 
with the SMA 2 JHo et alJEool on 28 April 2006. The data 
from both sidebands were combined, giving an effective band- 
width of 3 GHz. The SMA was in the compact configuration 
and baselines ranged from 5 to 52 metres. The low elevation 
of the sources as seen from Mauna Kea and the relatively short 
integration time of ~ 30 minutes per source resulted in a rather 
elongated beam of ~ 9 x 2.5 arcsec (natural weighting). The 



2 The SubMillimeter Array is a joint project between the 
Smithsonian Astrophysical Observatory and the Academia Sinica 
Institute of Astronomy and Astrophysics and is funded by the 
Smithsonian Institution and the Academia Sinica. 



raw visibility data were calibrated and flagged with MIR, and 
the calibrated visibility data were analysed with MIRIAD. The 
gains were calibrated on the quasar PKS 1622-297, and the 
absolute fluxes and correlator passbands were calibrated on 
Uranus with an expected uncertainty of ~20 %. All eight anten- 
nae of the array were available at the time of the observations. 
For an overview of the SMA observations, see Table|2] 

3. Results 

3. 1. Source fluxes 

From the 15 sources that were observed with the ATCA in 2003 
and 2005, ten were detected at better than 3<x, with 3. 3 -mm 
fluxes ranging from ~ 6 mJy up to ~ 30 mJy. For those ob- 
jects, both point sources and circular Gaussians were fitted in 
the (m, v) plane, the results of which are presented in Table [3] 
Those sources that were detected at better than 5<x are shown 
in Fig.^ overplotted on 2MASS K-band (2.1 /mi) images. The 
positions of the infrared sources and the millimetre peaks agree 
very well. 

Table 0] presents the results of the SMA 1.4-mm observa- 
tions of the Lupus sources, which were all detected at better 
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Table 2. Overview of the observations. 



Obs. dates 


Freqs. covered {u, v) range covered 


# Antennae used 




(GHz) 


(U) 




ATCA 


10 July 2003 


89.999, 90.095 


8-33 


3 


11 July 2003 


89.999, 90.095 


9-33 


3 


12 July 2003 


89.999, 90.095 


6-33 


3 


13 July 2003 


89.999, 90.095 


6-33 


3 


19 August 2005 


89.181, 91.456 


19-76 


5 


24 August 2005 


89.180,91.456 


12-70 


5 




90.000, 90.096 


14-75 


5 


25 August 2005 


89.176,91.456 


10-71 


5 


26 August 2005 


89.176,91.456 


12-71 


5 


28 August 2005 


89.176,91.456 


12-71 


5 


SMA 


28 April 2006 


217.347, 226.892 


4-53 


8 



Table 3. Basic results of ATCA observations at 3.3 mm. 



Source name 


Continuum flux" 


RMS b 


Gaussian size 


RA a 


Dec." 


Beam size c 


Effective frequency 






(mJy) 


(mJy/beam) 


(arcsec) 


(J2000) 


(J2000) 


(arcsec) 


(GHz) 




(P) 


(G) 














SYCha 




< 4.8 d 


1.6 




10 56 30.4 


-77 11 45.0 


2.5 x 2.2 


90.3 


CR Cha 


6.2 


6.1 


1.5 


< 1 


10 59 06.9 


-77 01 39.7 


2.5 x 2.1 


90.3 


CS Cha 


5.9 


8.8 


1.5 


1.65 ±0.96 


1 1 02 24.9 


-77 33 35.9 


2.4 x 2.2 


90.3 


DI Cha 




< 3.9 d 


1.3 




11 07 21.6 


-77 38 12.0 


2.4 x 2.2 


90.3 


KG 28 




< 3.9 d 


1.3 




11 07 57.9 


-77 38 50.0 


2.4 x 2.3 


90.3 


Glass I 




< 3.0 d 


1.0 




11 08 15.1 


-77 33 59.0 


2.4 x 2.4 


90.3 


KG 49 


11.9 


29.5 


1.5 


2.78 ± 0.40 


11 08 38.6 


-77 43 52.1 


2.4 x 2.2 


90.3 


WW Cha 


25.9 


33.1 


1.2 


1.32 ±0.16 


11 10 00.0 


-76 34 58.0 


2.5 x 2.2 


91.0 


XX Cha 




< 4.5 d 


1.5 




1111 39.7 


-76 20 21.0 


2.4 x 2.2 


90.3 


TCha 


6.4 


7.0 


1.0 


0.78 ± 1.03 


11 57 13.6 


-79 21 31.7 


2.5 x 2.4 


90.3 


HT Lup 


8.3 


12.0 


1.1 


1.40 ±0.41 


15 45 12.9 


-34 17 30.8 


2.4 x 1.7 


90.0 


GWLup 


8.5 


9.5 


1.9 


0.98 ± 1.22 


15 46 44.7 


-34 30 36.0 


5.3 x 1.7 


90.0 


IMLup 


8.9 


13.1 


1.3 


1.40 ±0.39 


15 56 09.2 


-37 56 06.3 


2.3 x 1.7 


91.5 


RULup 


12.7 


15.6 


1.0 


0.99 ± 0.32 


15 56 42.3 


-37 49 16.0 


2.3 x 1.7 


90.0 


HK Lup 


7.3 


9.1 


2.1 


2.85 ± 1.70 


16 08 22.5 


-39 04 46.3 


ill defined e 


90.1 



NOTE. — Source name in boldface indicates that the source is resolved at 3.3 mm with the ATCA, see Sect. 13.21 
a Continuum flux and position are from fits in the (u, v) plane. For sources that were detected at 3cr, both the point-source flux (P) and the 

integrated flux for a Gaussian (G) are shown. For sources that were not detected, the coordinates of the phase centre are quoted. 
b Calculated from the cleaned image. 
c Restored beam, using natural weighting. 
d Quoted value is 3cr upper limit. 

e Large, elongated beam due to short integration time and availability of only three antennae with 3-mm receivers. 



than 15<x, with fluxes ranging from ~ 70 mJy up to ~ 210 mJy. 
For comparison, the SE ST 1.3 mm single-dish fluxes from 
Niirnberge r et alJ ll 19971) are also shown. 

3.2. Are the sources resolved? 

A plot of the visibility amplitude as a function of baseline for 
each of our target sources is presented in Figs. [2] and [3] An 
amplitude that decreases as a function of baseline indicates a 
resolved source, suggesting that at least some of our sources 
are resolved. Note that in principle the resolved structure could 
be the protoplanetary disk, or the remnant centrally-condensed 
envelope around it. One can distinguish between the two by 



looking in more detail at the shape of the amplitude vs. (u, v) 
distance curves. A disk shows a shallower profile at shorter 
baselines, levelling off to the flux of the integrated disk emis- 
sion. A power-law envelope shows a flux that increases steeply 
towards the shortest baselines, and the total integrated flux can 
only be obtained with single-dish telescopes, observing with 
a larger beam. An illustrative example is provided by Fig. 3 
in Jeirgen sen et alJ ( 120051) . who need a disk and an envelope to 
explain the amplitude as a function of baseline for the deeply 
embedded class source NGC 1333 IRAS 2A. 

The amplitude vs. (u, v) distance plot for KG 49 (Cha IRN) 
appears to indicate a resolv ed envelope rather tha n a disk. This 
is in line with the results of Hennin g et al J 1 1 993 ) who find that 
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Fig. 1. ATCA images of the A = 3.3-mm continuum emission (contours), overplotted on 2MASS K-band (2.1 /jm) images 
(grayscales). Contour levels are drawn at 2, 4, 6, 10, 15, and 20 times the noise level; negative contours are dashed. The positional 
offsets are with respect to the fitted coordinates (see Table[3}. 



the SED from 1 pm to 1 mm cannot be fitted with a disk, but 
is consistent with a spherically symmetric model with constant 
density and Ay = 25 mag. Furthermore, all other sources show 
silicates in emission around 10 pm, which is indicative of a disk 
without significant foreground absorption. KG 49 on the other 
hand does not show silicate emission around 10 pm. Ho wever, 
it does show bands of H?Q, C O, and CO2 ice JGiirtler et alJ 
1 19991 Pontonnid an et alJ l2003). which can be explained by a 
cold envelope. It thus appears that KG 49 is not a genuine T 
Tauri star, but a less evolved object that is still embedded in 
an envelope. We will therefore disregard KG 49 in the further 
analysis in this work. 

We consider a source to be resolved if the integrated flux 
of the fitted Gaussian is at least 2cr higher than the flux ob- 
tained from a point-source fit (see Tables [5] and |4j. Note that 
for a point source the peak flux density equals the integrated 
flux. According to this definition, five of the detected sources 
(KG 49, WW Cha, HT Lup, IM Lup, RU Lup) are resolved by 
the ATCA, and two of the Lupus sources (IM Lup, HK Lup) are 
resolved by the SMA. In the case of KG 49 it is the envelope 
that is resolved, in the other cases the disks. 



3.3. Opacity index 

At low frequencies, i.e., in the Rayleigh-Jeans regime, the flux 
density, F v , is related to frequency, v, by a power law: F v oc v° . 
One can determine the dust opacity index, ft, where k v oc 
from the observed spectral index, a, through 



0*(or-2)(l + A), 



(1) 



where A is the ratio of optically thick to opti cally thin emission 
from the disk jBeckwith et al.ll990llBeckwith & Sargentll99ll 
iRodmann et alfcOOd) . 

Due to the frequency dependence of the dust opacity, pro- 
toplanetary disks are generally optically thick at short wave- 
lengths and optically thin at long wavelengths. Optically thick 
emission cannot be neglected in the inner disk where column 
densities get very high, even at very long wavelengths. The ra- 
tio of optically thick to optically thin emission coming from the 
disk is given by 



A*-/>x[(2-«)ln{(l-p/2)f}]- 



(2) 



where p and q are the disk's surface density and temperature in- 
dices, and f is the average disk opacity at the frequency under 
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Table 4. Basic results of SMA observations at 1.4 mm. For comparison, the SEST 1.3 mm single-dish fluxes from 
iNiirnberger et alJ i 19971) are shown. 



Source name 


Continuum flux 3 


RMS" 


Gaussian size 


RA a 


Dec. 3 


Effective frequency 


SEST 1.3 mm c 






(mJy) 


(mJy/beam) 


(arcsec) 


(J2000) 


(J2000) 


(GHz) 


Flux 


RMS 




(P) 


(G) 












(mJy) 


(mJy) 


HTLup 


73 


77 


4.0 


1.01 ±0.66 


15 45 12.9 


-34 17 30.1 


221.3 


135 


15 


GWLup 


64 


70 


3.7 


1.22 ±0.48 


15 46 44.8 


-34 30 35.7 


221.3 


106 


18 


IMLup 


188 


214 


4.3 


1.33 ±0.20 


15 56 09.2 


-37 56 06.5 


221.3 


260 


9 


RULup 


148 


158 


4.5 


1.02 ±0.32 


15 56 42.3 


-37 49 15.9 


221.3 


197 


7 


HK Lup 


89 


101 


3.9 


1.43 ± 0.38 


16 08 22.5 


-39 04 47.5 


221.3 


84 


17 



NOTE. — Source name in boldface indicates that the source is resolved at 1.4 mm with the SMA, see Sect. 13.21 
a Continuum flux and position are taken from fits in the (u, v) plane. Both the point-source flux (P) and the integrated flux for a Gaussian (G) 

are shown. 
b Calculated from t he cle aned image. 

c iNurnberger et alJ<1997l) . The SEST fluxes are in general higher than the SMA fluxes, partly due to the slightly shorter effective wavelength, 
1.3 mm (SEST) vs. 1.4 mm (SMA). 



Table 5. SEST 1.3 mm fluxes from the literature, mm slopes, 
derived opacity indices, and disk masses. 



Source 


SEST 1.3 mm 11 


Mm 


Opacity 


Disk mass c 


name 


Flux 


slope b 


index b 


A/ d , s k 




(mJy) 


a 


P 


CMo) 


SY Cha 


< 172 






< 0.011 


CR Cha 


125 ± 24 


3.2 ±0.5 


1.5 ±0.6 


0.014 


CS Cha 


128 ± 46 


2.9 ± 0.5 


1.0 ±0.6 


0.021 


DI Cha 


38 ± 11 


> 2.4 


> 0.5 


< 0.009 


KG 28 








< 0.009 


Glass I 


70 ±22 


> 3.4 


> 1.7 


< 0.007 


WW Cha 


408 ± 29 


2.7 ±0.7 


0.8 ±0.8 


0.077 


XX Cha 


< 252 






< 0.011 


TCha 


105 ± 18 


2.9 ±0.5 


1.1 ±0.6 


0.003 


HTLup 


135 ± 15 


2.5 ± 0.4 


0.4 ±0.5 


0.025 


GWLup 


106 ± 18 


2.4 ± 0.4 


0.5 ±0.5 


0.019 


IMLup 


260 ±9 


3.2 ± 0.5 


1.4 ±0.5 


0.027 


RU Lup c 


197 ±7 


2.5 ±0.1 


0.5 ±0.1 


0.032 


HK Lup 


84 ± 17 


2.5 ± 0.4 


0.7 ±0.5 


0.033 



iHenning et"al]h993h : fNurnberger et alj<1997l) . 
Calculated from SMA and ATCA fluxes (Gaussian fits) from this 
work and SEST fluxes, where available. 

Disk masses estimated from ATCA fluxes (Gaussian fits), assum- 
ing a gas-to-dust ratio *P = 100, a dust opacity at 3.3 mm k v = 0.9 
cm 2 g -1 , and a dust temperature T^, = 25 K. 
Includes JCMT fluxes from Weintrau betafl ll989l) . 



consideration (seelBeckwith et alJl990l for details). Following 
iBeckwith etafl Jl990l) and iRodmann etaD fcOOd) . we used 
p = 1.5. q is uniquely determined from the spectral index at 
wavelengths where the dust opacity is high {A < lOO/vm), and 
from IRAS photometry we obtained values ranging from 0.4 to 
0.7. f can in princi ple only be d etermined if the physical disk 
radius is known. Rodm ann et alJ 12006) use f = 0.01 at 7 mm. 
Taking j3 — 1 as a fiducial value for the opacity index, we adopt 
f = 0.02 at 3.3 mm and find values for A ranging from 0. 18 to 
0.22. 



The spectral index a and the opacity index /3 were deter- 
mined from the ATCA 3.3-mm fluxes and SEST 1 ,3-mm fluxes 
from the literature (tNiirnberger et al.l Il997l IHenning et alJ 
1993). The uncertainties in a and /3 are typically 0.5, due to 
both the large uncertainties in the absolute fluxes of the data 
points and to the relatively short wavelength range over which 
a was determined. For the Lupus sources, a was also deter- 
mined with the SMA 1.4-mm fluxes included, giving consis- 
tent results. The robustness of the value for a is illustrated 
by Fig. 0] which shows the integrated fluxes from this work 



and single-dish (sub)millimetre 


fluxes from the literature for 


RU Lup t Weintraub et alJ 


1989; 


Niirnbergeret al. 1997), along 



with the fitted slope: a = 2.5 +0.1. We thus find that the formal 
uncertainties in a and /3 for the other sources are probably over- 
estimates. Unfortunately, RU Lup is the only source for which 
currently observations at such short wavelengths are available. 
The low error in the fit to its data point illustrates that more 
observations, preferably over a larger wavelength regime, will 
help to better constrain values for a, and hence (3. 



The values for a and /? for the sources in our sample are 
presented in Table [5] Fig. |5] shows the cumulative fraction of 
source s with given dust-opacity indices for our sample and 
that of Rodma nn et alJ ([2006 ) . Rodma nn et alJ corrected for the 
contribution of free-free radiation at 7 mm. We estimate the 
contribution of free-free radiation to be ~ 5% at 3.3 mm, and 
hence it is not necessary to correct for it at this wavelength. 
Fig-EiUustrates that contamination from free-free emission is 
not an issue for RU Lup. The Kolmogoro v-Smirnov test gives 
a probability of 58% that Rodma nn et all s and our groups are 
drawn from the same distribution. However, the uncertainties 
in [5 are still quite large, and are not taken into account by the 
standard Kolmogorov-Smirnov test. 
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Fig. 4. F v vs. v for RU Lup. Filled symbols s how fluxes from 
single-dish observations (JC MT and SEST, IWeintraub et alJ 
1989; Niirnbe rger et al]ll997l) . open symbols show the values 
for the interferometric observations from this work. All points 
fall on a line with a = 2.5 + 0.1, where F v oc y° '. This suggests 
that there is no extended emission over the size scale of the pri- 
mary beams of the different telescopes, as may be the result of, 
e.g., a remnant envelope or ambient cloud material. Hence, no 
extended emission is filtered out in the interferometric observa- 
tions. 
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function. We used the integrated fluxes from the Gaussian fits, 
and assume d ¥ _ 100, T^ t — 25 K, and k y = 0.9 cm 2 g -1 at 
3.3 mm rcf. lReckwithetalJ ll990) 3 ]. The derived disk masses 
for our sources range from < 0.01 up to 0.08 M , and are pre- 
sented in Table|5] Note that the mass estimates are quite uncer- 
tain and may easily be off by a factor of 2, due to the uncer- 
tainties in the parameter values. Given the uncertainties in the 
measured quantities it is not meaningful to extend the analy- 
sis to include the contribution of optically thick emission and 
detailed disk structure. 

3.5. Molecular line emission 

Both IM Lup and WW Cha were observed in spectral-line 
mode to search for HCO + J = 1-0 emission. In general, CO is a 
more easily detectable gas tracer in the millimetre regime, but 
since the ATCA does not presently have the capability to ob- 
serve at frequencies above 106 GHz, HCO + was used to inves- 
tigate the gas component of disks. We tentatively detect HCO + 
at ~ 0.05 Jy beam -1 in IM Lup, as presented in Fig. [6] Van 
Kempen et al. (2006) detected a double-peaked feature, consis- 
tent with a rotating gas disk, using JCMT observations of the 
12 CO J = 3-2 line. However, they find the line at ~1 km s _I 
lower velocity. The cause of this discrepan cy is unclear. 

A ray-tracing programme (Hogerheiid e & van der Takl 

2000) was used to compute the line profile of HC O + 7=1- 
0, usin g the model for IM Lup described in Ivan Kempen et all 
ll2006h . The abundance of HCO + was set to be 10" 8 with re- 
spect to H2 in areas with temperatures above 30 K, and 10~ 12 in 
areas with temperatures below 30 K. An intensity of ~ 0.02 Jy 
beam -1 was predicted, consistent with the observations within 
the uncertainties. 

No HCO + line emission is detected to a limit of 97 mJy 
beam 1 (3<x) in the d irection of WW Cha . Combining our new 
results with those of W ilner et alJ 1120031) who detected HCO + 
in TW Hya, but not in HD 100546, we have two detections 
of HCO + and two non-detections. Data for more targets are 
required, in order to determine whether the presence of HCO + 
is related to the disk's evolutionary state. 



Fig. 5. Cumulative number of sources with an opacity index less 
than a given val ue of B for the sources i n our sample (solid line), 
those studied by Rodmann et al. (2006) (dashed line), and those 
studied bv lNatta et alJ j2004) (dotted line). 



3.4. Disk masses 

For optically thin disks, the dis k mass is directly pro- 

porti anal to the flux F v (see, e.g., Hildebra ndll983l[Natta et all 



M disk = 



F v y>D 2 
K V By(T dusl y 



4. Discussion and interpretations 

4. 1. Grain growth 

About half the sources that are detected are spatially resolved. 
The inferred physical sizes (disk radii of ~ 100 AU) indicate 
that the emission coming from the disks is predominantly op- 
tically thin at millim etre w avelengths. This is illustrated by the 
work of Testi et al. (2001), who do detailed modeling of the 
Herbig Ae stars UX Ori and CQ Tau, and find that disk radii 
of ~ 100 AU, combined with millimetre slopes of < 3 and 
1.3 mm fluxes of ~ 10 2 mJy, are well explained by an opac- 
ity index B = 1 .0. If we assume that all sources in our sample 

^ 3 Note that Beckwith et al. 1 199G) estimate k v = 0.1(v/10 12 Hzf 



where T* is the gas-to-dust ratio, D is the distance to the source, 
k v is the dust opacity, and B y (T,i usl ) is the brightness at the fre- 
quency v for a dust temperature T dust , as given by the Planck 



cm 2 g _I , where k v is the opacity index for the gas and the dust com- 
bined, i.e., with an implied gas-to-dust ratio. Our values for k v , how- 
ever, are for the dust alone, and hence we have to account for the 
gas-to-dust ratio explicitly. 
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Fig. 6. Spectrum of HCO + J = 1-0 line emission observed from 
IM Lup, at the position of the peak continuum emission, binned «? 
to 0.44 km s velocity resolution, in the 2'.'3 x I'.'l synthesized 
beam. » 1-0 



are optically thin at millimetre wavelengths, the opacity index 
P can be determined. Seven out of ten detected sources have 
P < 1, which can be naturally explained by grain growth to 
millimetre and centimetre sizes (Draine 2006). 

4.2. Comparison with Spitzer infrared data 

Our values for ft indicate grain growth to sizes of millimetres 
and larger in the outer disks. The surface layers of the inner 
disks can be probed by infrared observations and, as noted 
in Sect. [2 the 10-yum feature indicates the growth of grains 
through changes in its strength and shape. 

In Fig. we compare the millimetre slope a as derived in 
Sect.[5]with the "strength" or peak 10-jum flux (Sj^ak) and to 
the "shape" or the ratio of the 11.3 to 9.8 //m flux (Sn.3/Sg,s) 
of the 10-yum silicate feature for those sources in our sa mple 

^^^^^^^^^^^^^^^^^^^^^^^^^^ ^^^^^^^^1 

that overlap with the samp les of Przvgod da et alJ y003) and 
iKessler-Silacci et aflfcOOcjl) . Here the normalized 10-//m spec- 
tra S v are given by 

(Fy ~ Fyj.) 



1 + 



< Fy, C > 



(4) 



where F Y is the observed spectrum and F Y>C is the fitted contin- 
uum, and < F y>c > is the frequency-averaged continuum flux 
(see Kessler-Silacci et al. 2006 for details). We find a positive 
correlation between a and Slj^ and a negative correlation be- 
tween a and S 11.3/59,8, especially if we leave out the source 
T Cha, which shows emission from polycyclic aromatic hydro- 
carbons at 11.3 fim and may be considerably older than the 
other sources in the sample. As grains grow from submicron 
sizes to several microns, the 10-Atm feature become s weaker 
and less peaked (see, e.g jKessler-Silacci et al.l2 006), and thus 
the correlations we find may well indicate that once grain 
growth starts, the grains quickly grow from (sub)micron sizes 

to millimetre sizes and larger, both in the inne r and outer disks. 

1 I 

This is consistent with grain-growth models ( Weidenschilling 
1988l ll997HDullemond & Dominikl2o"o4bll2005l) . which show 
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Fig. 7. Slope in the millimetre regime vs. the peak 10-yum flux 
(upper panel) and the ratio of the 11.3 to 9.8 /mi flux (lower 
panel) of the 10-/im silicate feature, for those sources in our 
sample for which data are available. 1=HT Lup, 2=GW Lup, 
3=IM Lup, 4=RU Lup, 5=CR Cha, 6=WW Cha, 7=Glass I. The 
source T Cha is depicted with a diamond. 



that growth to metre sizes is rapid (~ 10 3 yr at 1 AU; ~ 10 5 yr 
at 30 AU). 

On the other hand, while a significant fraction of the dust 
has already grown to sizes of millimetres and centimetres, the 
infrared data indicate that micron-sized grains are also still 
present, at least in the surface layers. This may be explained 
by assuming that not only grain gro wth, but also aggregate 
fragm entation takes place in the disks dDullemond & Dominild 
l2005h . After about 10 4 years a semi-stationary state is reached 
for sizes below a < 1 cm, which may last for several 10 6 years. 

Note that Fig.|7]shows a trend rather than a bimodal distri- 

1 n il - ■ ■ 

bution. Dullemond & Dominik (2005) suggest fragmentation 
of grains to allow for the semi-stationary state that is observed 
in the disks around T Tauri stars. If the correlation between the 
10-yum feature and the millimetre slope is confirmed by a more 
extensive dataset, it would indicate that when aggregates are 
fragmented in collisions, the size of the fragments increases as 
the aggregate sizes increase. Hence, when the largest particles 
grow from millimetre to centimetre sizes, the submicron-sized 
grains are no longer replenished, and the 10-/im feature - trac- 
ing the upper layers of the disk where the stirred-up small par- 
ticles reside - flattens while the millimetre slope - which traces 
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the midplane where the largest particles are present - becomes 
shallower. 

We also compared a to the spectral index o ver the 13 
to 35 /mi range from iKessler-Silacci et.alJ J2006). However, 
a could only be d etermined for five of the sources in 
IKessler-Silacci et all and no obvious correlation is found for 
this small sample. 

4.3. Comparison with Herbig Ae/Be stars 

It is interesting to note that lAcke & van den Anckerl d2004l) do 
not find a correlation between the 10-/mi silicate feature and the 
(sub)millimetre spectral index f or their sample of 26 Herbig 
Ae/Be stars. On the other hand, Ack e et al.l (|2004) do find a 
correlation between the shape of the mid-IR (12-60 fim) SED 
and the (sub)millimetre slope in their sample. They interpret 
this finding as a correlation between disk geometry and grain 
size: as grains grow, the dis k structure evolve s from flared to 
geometrically flat (see, e.g.. lDullemondlE002l) . For the twelve 
sources in our sample for which it was possible to determine a, 
we determined the non-colour-corrected//M5 [12]-[60] colour 
and compared this to a. No clear correlation was found, con- 
sistent with the above lack of correlation with the 13-35 fim 
spectral index. This may indicate that grain growth has a less 
dramatic effect on the geometry of disks around the low-mass 
T Tauri stars than on the intermediate-mass Herbig Ae/Be stars, 
or that IRAS fluxes are contaminated by other (extended) emis- 
sion in the larg e IRAS beams. 

iNatta et alJ {2004) did a study similar to ours for a sample 
of six isolated Herbig Ae stars and three T Tauri stars. They 
estimated /3 in two different ways, firstly using only interfero- 
metric data, and secondly with single-dish data included. Their 
results are included in Fig.|5] where the values for/? with the 
single-dish data included were used for comparison. According 
to the Kolmogoro v-Smirnov test there is a 65% probability 
that the sample of Nat ta et all and our sample are drawn from 
the same distributi on. It is interesting to note that in general 
the sources in Natt a et all s sample are conside rably older than 
those sources in the current work and that of Rodmann et al. 

This may be an extra indication that the size distri- 
bution of grains in protoplanetary disks remains stationary for 
several Myr once sizes of millimetres are reached. 

5. Concluding remarks 

We used the ATCA to make interferometric observations of fif- 
teen southern T Tauri sources at 3.3 mm. Ten of the sources 
are located in Chamaeleon, the remaining five in Lupus. The 
Lupus sources were also observed at 1.4 mm with the SMA. 
The main results are as follows. 

- Five of the Chamaeleon sources were detected at the 3<x 
level or better with the ATCA, for the other five sources we 
have strict upper limits. 

- All five Lupus sources were detected at better than 3cr with 
the ATCA, and at better than 15<x with the SMA. 

- Five of the ten detected sources are spatially resolved with 
the ATCA, and two of the five Lupus sources are resolved 



with the SMA. This indicates that the emission coming 
from these sources is optically thin. 

- Adopting optically thin emission for all sources in our sam- 
ple, we estimated the opacity index /?. We find that /3 < 1 
for seven of the ten detected sources. Such values for ft 
can be naturally explained by grain growth up to millimetre 
sizes and beyond. 

- We find a tentative correlation between the millimetre slope 
and the peak IQ-fim flux and ratio of the 11.3 to 9.8 fim 
flux of the 10-fim silicate feature. This indicates that grain 
growth takes place in the outer disk and in the surface lay- 
ers of the inner disk simultaneously. It also confirms earlier 
findings that agglomerate destruction must take place in the 
disk to preserve the small grain population. If confirmed by 
more extensive datasets, this correlation may be an indi- 
cation that larger aggregates produce larger particles when 
fragmented in collisions. 

- HCO + was tentatively detected in IM Lup, and was not seen 
in WW Cha. 

This type of work will be greatly assisted by the upcoming 
7-mm upgrade for the ATCA, as well as by the installation 
of a new correlator via the Compact Array broadband back- 
end project. Estimating the millimetre slope over a wavelength 
range up to 7 mm will reduce the error in a - and hence in [5 
- from ~ 0.5 to ~ 0.2. At these longer wavelengths, however, 
the contribution of free-free emission becomes important, and 
lower-frequency measurements are required to separate the var- 
ious contributions. A pilot study to determine the contribution 
of free-free emission at centimetre wavelengths in our sources 
was carried out 

The capability to observe the continuum across a 2 GHz 
bandwidth, as opposed to the current 128 MHz, will provide a 
marked increase in sensitivity, allowing more disks to be de- 
tected. Furthermore, as the new generation of (sub)millimetre 
interferometers, notably eSMA, CARMA, and ALMA, be- 
comes available to the scientific community, it becomes pos- 
sible to resolve the circumstellar disks at millimetre and sub- 
millimetre wavelengths down to subarcsecond scales. A sur- 
vey of a large sample of objects, in different clouds, and at 
higher resolution will greatly enhance our understanding of the 
timescale for grain growth and the building of planetary sys- 
tems as a function of disk radius. 
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Fig. 2. Amplitude vs. (u, v) distance for sources observed with the ATCA. The data points give the vector-averaged amplitude per 
bin, where the data are binned in annuli according to (u, v) distance. The error bars show the statistical lcr errors, and the dotted 
lines give the expected amplitude for zero signal. 
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3. Same as Fig.|2]for sources observed with the SMA. 



